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The standard cosmological model, which includes only Compton scattering photon interactions 
at energy scales near recombination, results in zero primordial circular polarization of the cosmic 
microwave background. In this paper we consider a particular renormalizable and gauge-invariant 
' standard model extension coupling photons to an external vector field via a Chern-Simons term, 

which arises as a radiative correction if gravitational torsion couples to fermions. We compute the 
transport equations for polarized photons from a Boltzmann-like equation, showing that such a 
coupling will source circular polarization of the microwave background. For the particular coupling 
\ considered here, the circular polarization effect is always negligible compared to the rotation of 

. the linear polarization orientation, also derived using the same formalism. We note the possibility 

that limits on microwave background circular polarization may probe other photon interactions and 
related fundamental effects such as violations of Lorentz invariance. 
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I. INTRODUCTION 

One of the great successes of the standard cosmology is the prediction and measurement of the temperature 



anisotropies in the cosmic microwave background radiation. Most of these photons have freely propagated since 
the epoch of last scattering roughly 14 billion years ago and encode the initial conditions for structure formation. 
Measurements are now consistent to high precision with the simplest cosmological models with an initial power-law 
spectrum of adiabatic perturbations. Linear polarization of the microwave background fluctuations is also a generic 
result of these models; recent detections of the linear polarization power spectrum and of the cross-correlation between 
linear polarization and temperature are also consistent with the same cosmological models. 
, In general, radiation can have linear polarization, with two degrees of freedom (a polarization amplitude and 
■ orientation) as well as circular polarization, with a single degree of freedom. It is well known that if an initially 
unpolarized photon field evolves solely via Compton scattering from free electrons plus free streaming, the resulting 
radiation field can have linear but not circular polarization. In the tight-coupling regime prior to last scattering 
when Compton scattering is rapid compared to the cosmological expansion time scale, the cosmic radiation field 
will be unpolarized. As the universe cools and the free electrons become bound into neutral hydrogen, a small 
linear polarization is generated from the balance of free-streaming and Compton scattering during this recombination 
. \ process, but the resulting microwave background radiation today has circular polarization which is identically zero. 
In this paper, we consider a generic class of interactions between photons and an external field which can produce 
circular polarization. The interactions have been considered in other contexts and are general enough to be expected 
in broad classes of theories beyond the standard model of particle physics. The same interaction can also arise if 
non-zero spacetime torsion impacts the microwave background radiation. The goal of this paper is two-fold. First, 
we provide an explicit calculation showing how circular polarization can be generically sourced in the microwave 
background, with the relevant evolution equations. Second, we demonstrate what the underlying microphysics might 
look like. 

Consider the following extension of the photon sector of quantum electrodynamics: 

•£ = ^MAXWELL C.T 

= —F^F^ + ge^A^F^ (1) 

where Ct is CPT odd and violates Lorentz invariance and g is the coupling constant of the interaction. Several authors 
have investigated such a Lorentz-invariance violating extension of QED for a constant 4- vector (see e.g. [H, [3,0] 
and references therein). These so-called Standard Model Extensions have been shown to be renormalizable while 
maintaining gauge invariance We consider here only the flat-spacetime interaction term Ct for simplicity; this 
will be a good approximation in any cosmological context (see Refs. J|[ and @ for the curved spacetime generalization, 
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which includes an extra factor of the square root of the metric determinant). If T M is fixed as a constant, the Lagrangian 
density in Eq. ([T]) is U(l) gauge-invariant apart from a boundary term; therefore in these cases the action is gauge 
invariant. 

It is well known that such an extension should result in optical activity in the propagation of electromagnetic 
radiation 0, H, [E E3] ; specifically, a modification to the dispersion relations of free electromagnetic radiation results 
in a rotation of the plane of linear polarization during propagation. We make no explicit assumption about whether 
is spacelike or timelike, although the timelike case appears pathological since it leads to a violation of causality 
and unitarity . The magnitude of optical activity of electromagnetic radiation has been constrained by anal ysis 
of observational data from cosmological sources and from the microwave background radiation 0, EE EE El, EE HE 
EE El) EE H3| ■ Such a term may arise as a radiative correction following the coupling of gravitational torsion with 
fermionic matter [U • The same term has also been associated with the cancellation of gauge anomalies in QED 
when the background field T p is allowed to couple to the axial current (see e.g. [HI). 

In this paper, we show that in addition to the well-known polarization rotation, such a term may also generate 
circular polarization, although for the specific case of Eq. ([1]) the circular polarization is always negligible compared 
to the polarization rotation. The generation of circular polarization following the optical activity produced by 
parallels the Faraday conversion and Faraday rotation effects for propagation in magnetized plasmas; for a discussion 
in the context of the microwave background, see [2~ij . The observation of circularly polarized microwave background 
radiation could be evidence of Lorentz-invariance violation and thus physics beyond the standard model; conversely, 
limits on circular polarization may constrain a certain class of standard model extensions. For a related analysis 
using an axion-like pseudoscalar coupling to the electromagnetic field, see Ref. (25j . who also find a nonzero circular 
polarization and rotation of linear polarization. 

In Sec. UH we review the usual description of polarized electromagnetic radiation in terms of Stokes parameters; 
linear polarization is described by the Q and U parameters, while circular polarization is described by a non-zero 
V parameter. Section IIIII reviews the construction of the Boltzmann-type equation for the photon number density, 
starting from the quantum-mechanical evolution of the photon density matrix. In Sees. IIVI and [Vj we calculate the 
evolution of the Stokes parameters to first and second order, respectively, in the interaction term Ct, deriving the 
evolution equation for the V polarization, which is generated from linear polarization due to the interaction term. 
We conclude in Sec. I VII with estimates of the size of the V polarization in the microwave background for given 
interactions along with the magnitude of linear polarization rotation. The mathematical details of evaluating the first 
and second-order interaction terms are relegated to Appendices A and B, while Appendix C addresses the issue of 
gauge invariance. 



II. STOKES PARAMETERS 



The polarization state of light is most easily described by making use of the Stokes parameters. For a complete 
review see, e.g., Refs. [IE HE HH or any optics text. Here we review the basic construction of the Stokes parameters 
in the classical and quantum mechanical contexts in order to motivate the quantum field theoretical construction. 
Consider a classical electromagnetic plane wave with electric field given by the components 

E\{£) = ax sm.(u)t — ei) and ^(t) = 0,% sin(a;i — £2) (2) 

where we assume, for simplicity, that the wave is nearly monochromatic with frequency u>, such that a±, 02, ei, and 
£2 only vary on time scales long compared to The Stokes parameters in the linear polarization basis are then 

defined as 

I EE (Kf + M 2 ), (3) 

Q ee (( fll ) 2 -M 2 ), (4) 
U ee (2aid2 cos 8) , (5) 
V ee (2aid2 sin 8) , (6) 

where 8 = £2 — £1 and the brackets signify a time average over a time long compared to uo^ 1 . The / parameter measures 
the intensity of the radiation, while the parameters Q, U , and V each carry information about the polarization of the 
radiation. Unpolarized radiation is described by Q = U = V = 0. The linear polarization of the radiation is encoded 
in Q and U, while the parameter V is a measure of elliptical polarization with the special case of circular polarization 
ocurring when a\ = ai and 8 = ±7r/2. From here on we will simply refer to V as the measure of circular polarization, 
which is technically correct if Q = 0. Note that while / and V are coordinate independent, Q and U depend on the 
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orientation of the coordinate system used on the plane orthogonal to the direction of propagation. Under a rotation 
of the coordinate system by an angle <f), the parameters Q and U transform according to 

Q' = Q cos(20) + U sin(20), 
U' = -Q sin(2</») + U cos(2<£), 



while the angle defined by 



~ 1 (U 
<s = — arctan I — 



goes to $ — <fi following a rotation by the angle (j>. Therefore, Q and U only define an orientation and not a particular 
direction in the plane: after a rotation by 7r they are left unchanged. Physically, this is simply a manifestation of the 
oscillatory behavior of the electric field. These properties indicate that Q and U are part of a second-rank symmetric 
trace- free tensor P a b, i.e. a spin-2 field in the plane orthogonal to the direction of propagation. Such a tensor can be 
represented as 



P — 



P 

-P 



in an orthonormal eigenbasis, where P = (Q 2 + J7 2 ) 1 / 2 is often called the magnitude of linear polarization. 
In quantum mechanics we can express the state of a photon A as 



(7) 



(8) 



where |ej) (i — 1,2) span the polarization state space and a% are in general complex. The projection operators 

I = ^(exl + h)^!, (9) 

Q = | ei )(e 1 |-|e 2 )(e 2 | ! (10) 

U | ei ) <e 2 | + |e 2 ) <ei| , (11) 

V = i\e 2 )(e 1 \-i\e 1 ){e 2 \ (12) 

have expectation values in single photon states which give the classical Stokes parameters, Eqs. ©-(JS]). In a general 
mixed state, the density matrix p on the polarization state space encodes the intensity and polarization of the photon 
ensemble. For example, 



(Q) 



1 



trjpQ) 

tr(p) tr(p) ' 



-tr 



Pll P12 
P21 P22 



1 

-1 



Pll - P22 

tr(p) 



Similar relations hold for the other "Stokes operators" such that the density matrix can be represented as 

MP) ( 1 + Q U-iV 
U + iV 1-Q 



(13) 



where Q = (Q), U = (U), and V = (V). 



III. THE PHOTON BOLTZMANN EQUATION 



We now review the construction of the evolution equation for the photon number operator under the influence of 
someperturbation to a free theory. The following formalism was developed to study neutrino mixing and damping 
[2^ . l30j ] . It has also been applied to describe the generation of linear polarization in the microwave background due 
to Compton scattering during recombination [3lj ] . generalizing an earlier kinetic equation treatment of microwave 
background temperature fluctuations [32j ■ 

Consider an ensemble of free photons. We will assume that the interaction Ct is slowly "turned on" and that the 
interactions of the photons with the external field T M are localized such that the photons can be considered free (with 
respect to the interaction Ct) both before and after each point interaction. That is, we make the usual assumptions 
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of scattering theory. We will not consider any possible interference effects which might occur between Ct and any 
other interaction. 

The free photon field in the Coulomb (radiation) gauge can be expressed as 

^3^0 [a s (k)e s „(k)e- ik -* + a\(k)e^(k)e^} , (14) 

where e sfjL {k) are the photon polarization 4-vectors and s, which takes the values 1 and 2, indexes the orthogonal 
transverse polarizations. The free creation and annihilation operators satisfy the canonical commutation relation 

Oi(fc), fit (A')] = (27r) 3 2fc°,5y(5 (3) (k- k') (15) 

where k° = |k|. 

We will be interested in the evolution of the polarization state of a photon ensemble, which is completely charac- 
terized by the density matrix pij defined via 

(fil(fc)fij(fe')) = (27r) 3 2fc°J (3) (k - k>y(k). (16) 

The number operator X\,(k) = fi](k)fij(k), according to (I16|) . is related to py(k) as 

(^■(k)} = (27r) 3 2fc° ( 5( 3 )(0K(k). (17) 

The infinite factor S^{0) is a remnant of the infinite quantization volume. As we show explicitly below, it cancels 
from all final expressions. Motivated by the construction in the quantum mechanical system above, we can project 
out quantities analogous to the classical Stokes parameters: 

tr(a p(k)) -> 7oc(W (k)> (18) 

tr(<7,p(k)) -> Qcx(P 11 (k))-(P 22 (k)) (19) 

tr{a x p(k)) -> C/oc(Pi 2 (k)) + (P 21 (k)) (20) 

tr(<r„p(k)) -» iV cx (P 12 (k)) - (£> 21 (k)> (21) 

where ao is the 2x2 identity matrix, at are the Pauli matrices and the trace is over polarization indices. 7io(k) is 
the free energy density operator of the mode with wavenumber k. 

We ignore any correlations such as (ai(k)aj(k')} and (fif (fc)at(fc')) which might be generated via the interaction Ct- 
In essence we are assuming that the background field T M varies slowly enough in time so that physical two-photon 
states are neither created nor destroyed by the interaction Ct- If we define ujt as a characteristic energy scale of the 
background field and lj as the frequency of a particular mode associated with free oscillations of the creation and 
annihilation operators, then we are requiring that lot <C Suj/uj <C to [33l |. Note that Suj/ui is the order at which mixing 
occurs between hi and fit, and such a mixing will result in a variation ^(fij(fc)fij(fe')) ~ (5lo / Lj)(a\(k)a,j(k')) , which 
again we neglect. Naively, we then expect the following formalism to fail for low-frequency photon modes, although 
at precisely what scale the approximation breaks down depends on the characteristic scale of the background field T M . 

The evolution of the number operator X>y (k) can be computed using a perturbative expansion in the interaction 
strength. The following reviews the construction detailed in [3(| ■ Recall that the time evolution of any Heisenberg 
picture operator is given by 

|p«=i[ff,C«]. (22) 

If the full Hamiltonian can be split into a free and interacting part, H = Hf ICC + Hint, then (|22[) becomes. 

d 



dt 



Vij = % [fff roc , Vij] + i [ffi„t, Vij\. (23) 



A first order perturbative approximation for the evolution of the number operator is given by replacing all operators 
on the right hand side of equation (|23p by their free theory counterparts, e.g. 



-V tJ )^ l {[Hl t ,V a ]), (24) 
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where O corresponds to the operator O evaluated in terms of the operators of the free theory. The above assumes 
that [Hf ICC ,T>°] = 0. We will refer to the term on the right hand side of (j2~4")l as the refractive term, or the forward 
scattering term. To determine a second order perturbative approximation, we will use the fact that we can expand 
any operator to first order in interactions as 

i(t) ~ i°(t) + i f dt'[H? at (t - f), i°(t)], (25) 

Jo 

with the initial conditions f (0) = |°(0). The expansion Eq. (|25p can be verified by explicitly taking the time derivative 
of both sides and seeing that one recovers the Heisenberg equation to first order in interactions. We now expand 
[flint, T>] as in (|2"5|), insert the result into (f2"5)) . and upon the evaluation of all operators in terms of the free theory 
operators arrive at 



dt 3 



(t)ci([H° nt (t),V ])- / dt'({H? nt (t-t'),[H? nt (t),V } 



(26) 



The second term in (|26[) will be referred to as the damping term or the non-forward scattering term. In terms of p, 
the evolution equation reads 



(2 7 r) 3 < 5 3 (0)2q |p u (0,q) = i([£P t (Q),£°.(q)]) 



dt' ([H? nt (t%[HU0)M(<i)}]), 



(27) 



where, as mentioned above, all factors of 6 3 (0) will cancel from the final expressions. In going from (|26|) to (|27|) . we 
have assumed the time step t in Eq. (|2"6"|) is both small relative to the characteristic time scale of the evolution of 
p and large relative to the duration of a single interaction. This allows us to take t — > oo and set p(t) = p(0) [30| . 
We have then replaced the integral (J °° dt) with (i dt), the difference being a principle part integral which is a 
second order correction to the refractive term. Equation (|27[) can be viewed as a Generalized Boltzmann Equation for 
the phase space function p. In this approximation, we have a set of differential equations for the components pij at 
t = 0; if the interactions are "forgotten" between intermediate collisions (an assumption known as molecular chaos 
in the derivation of the standard Boltzmann equation), then the differential equations will be valid for all times over 
which the interaction is relevant. The Liouville terms on the left side will incorporate any effects which result from 
a departure of the spacetime metric from a flat metric, including any weak inhomogeneities due to the presence of 
gravitational perturbations about a homogeneous cosmology. The case H lnt (t) — J eta/^ A^ip as in full QED, where 
tp is a spinor field associated with the electron, recovers the radiative transfer equations of Chandrasekhar [3J] in the 
appropriate limits (see [3l[ for more details). 

As we will see below, this construction allows linearizing the right side of Eq. ([27]) in p. We can expand the photon 
density matrix about a uniform unpolarized distribution (ignoring any small inhomogeneities) as 



p t] (t,k,k) = P f(t,k) + p\;>(t,k,k), 



(28) 



where p^ — p§l and p y yj — P21 — 0- As a consistency check, the right side of Eq. (|27[) should vanish when evaluated 
in terms of p^ so that in for example an FRW background of zero spatial curvature with scale factor a(t) we have 

m 



,(0) _ JO) _ 



d_ (o) _ dj^_ _ d flpW _ 
dt Pl1 ~ dt a dk ~ u ' 



(29) 



the solution of which is Pxi(t, k) = p^ ; (fca), recovering the uniform redshift due to cosmological expansion. 

In order to make contact with the measurable Stokes parameters, we define the normalized brightness perturbations 



J°) 



31 



A/ = 


yqdpf}(q)- 
U dq - 


"(piU&A 


(30) 


Aq = 


-4 dq - 


~\p$- P £h 


(31) 


A v = 


\qdpii (q)~ 

.4 dq - 


~\p$+p£h 


(32) 


A v = - 


.fqdpuiq)]- 1 , ji) ji)n 

J4 dq 1 [Pl2 P2lh 


(33) 
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where q = ka is the comoving photon momentum and we have expanded the density matrix p in a linear polarization 
basis. 



IV. THE FIRST-ORDER INTERACTION TERM 



We now evaluate the right side of the Generalized Boltzmann Equation (|27|1 for an interaction Hamiltonian which 
is linear in the Hamiltonian density 

Ht = -ge^ : A^T v F a0 : , (34) 

where A^ for the free theory is given by Eq. (fT4|) and F^ v = 2d[^A u ^ is the free electromagnetic field strength operator. 
We will treat T M as a classical background field, the dynamics of which are not influenced by the electromagnetic 
interaction Eq. ([34"]) and are assumed to be externally prescribed. The symbol : • • • : denotes normal ordering of the 
enclosed operator products. 

Following a local U(l) gauge transformation SqA^ = d^X of A^, the resulting change in the Hamiltonian density 
Eq. (|3~4"|) is given by 

SgHt = -ge^id^TvF^ = ge^ a PX{d^T u] )F a0 - ge^d^XT u F a0 ), (35) 

where the square brackets denote anti-symmetrization and we have used d^F a ^ = 0. If d^T^ = then the interaction 
Hamiltonian is gauge invariant apart from a boundary term. We assume that the external field T^(x) does indeed 
satisfy this condition allowing us to maintain gauge invariance, which is detailed in the second-order calculation in 
Appendix O Furthermore, we assume that T^(x) is a pseudo- vector: under a parity transformation the external field 
transforms as 

~ x) = { vr}{i xt 1 = 1, 2, 3 (36) 

and under time reversal the external field transforms as 

We will find it useful to consider the Fourier transform of T M , 

T„(p) - J d 4 x T^(x)e- tp - x . (38) 
In momentum space the gauge invariance restrictions can be expressed by the condition 



Pfrf v] (p)=0. (39) 



In order to determine whether circular polarization can be sourced at some order of (|35[) . we will not need to impose 
any further restrictions on T^x) aside from those listed above. Specifically, our calculations do not assume that T^x) 
is either timelike or spacelike. As indicated by the investigation in [Tlj . violations of causality and unstable solutions 
may arise for the case of timelike T^ix), and the results here must be interpreted with care in this case. 

In Appendix [5] below we detail the calculation of the first-order interaction Hamiltonian which is linear in the 
Hamiltonian density Eq. (|34|) . as well as the refractive and damping terms of Eq. (|27|) due to first-order processes. 
Quoting the results of Appendix [Al the interaction Hamiltonian for first-order processes is given by 

/T (nW^o+po — |k+p|)t 
dpdk 2|k + p| ^(fc)a r (fc)e: M (k)(e^(|k| + |k + p|)+e^(2^+^))e^(k+p). (40) 

The refractive term of Eq. (|2T[) for this first order interaction Hamiltonian is given by 

^[^(Oj.A.vCq)]) = 4gq O (2n)H^Q)e s ^q)(S urPsv (q)~5 vs p ur (q))A^(q)e r0 (q), (41) 
where _4 M/3 is defined via 

(2tt) 3 5 3 (0)A^(q) = (2n) 3 S s (0)A^(q) = J ^-^(2 7 r) 3 <5 3 (-p)f,(p)e^(2 (7tt + p a ), (42) 
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and we have for convenience defined p = (Ag, p), Ag = |q + p| — |q|. 

As detailed in Appendix [A] the damping term of Eq. P7| due to first order processes is given by 



D„ 



dt' 



-(2ig) 2 e mK e s 



dp 



-f x (p)q p e aX P K f y (p)(2q + p) a e^e n(T e r0 



2 q + P 

x I d nr 5 us p mv (q) + 5 rn S vs p um (q) - S us S mv p rn (q + p) - S um 5 sv p nr (q + p) 



(43) 



where we have defined e rfi = e rfJ _(q) and e rfl = e rfl (q + p)- 

Once the density matrix is expanded as given in Eq. (f28|) . it is straightforward to see that the refractive term, Eg . (|4Tj) . 
vanishes when evaluated in terms of p(°)(|q|). For the damping term, Eq. (|43|) . we must perform an expansion of 
p (0) (|q + Pl) = P (0) (|k|), in the p integral of Eq. l[4"5|). about |q|: 



p( )(|k|) = p(°)(|q|) + O(^-(|q| 



(44) 



which is a suitable approximation as long as T(p°,p) has support solely over |p| <C |q|, where |q| is the energy of the 
scattering photons. Then to lowest order in p^ the damping term vanishes and we have 



d_ 

dt 



d\q\ 



(45) 



In terms of the Stokes brightness perturbations defined in Eqs. ([3T))) -- (f3"3")) . to first order in g the evolution of the 
polarization of the photon ensemble becomes 



d a 
dt^ 
d a 
dt Au 
d \ 
dt Av 



0, 

-ga(q)A u , 

ga(q)A Q , 

0, 



where we have defined the quantity 



a (q)^4e^(q)A^(q)e 20 (q). 



(46) 
(47) 
(48) 
(49) 

(50) 



For processes which are first order in the Hamiltonian density Eq. (f3"4"|) , according to Eq. (|4"5|) no circular polarization 
is generated to O(g) in our approximation. In fact, it is easy to see that to O(g) (the refractive term), Eqs. |46|) 



fl4"9")) reproduce the well-known effect of optical activity of the electromagnetic radiation, rotating the plane of linear 
polarization during propagation 0, @, @, 101 . This is a useful check of the calculations. 

The relevant linear combinations of the damping term Eq. (f4"3")l which source the polarization of the photon ensemble 
to 0(g 2 ) and due to first order processes are 



Ai + A^ = -{2igf 



dp 



2|q + p| 



T x {p)q p e° x ^T u {p)(2q+p) a 



x \ {hah/3 + hahp) 



2ei K ei At /^ 1 1 ) (q) + 2e2 K e 2p p 22 ' (q) + {^i^ + eiK^lPv (<l) + P^fa)] 



(1), 



2e V ei CT/ oi 1 1 ) (q + p) + 2l 2& l 2(J p 1 ^ {q + P) 



+ (hpha + £2/3eia)[Pv2 (<l + P) + P21' + P)] 



0\g 



2 V 0) (|q|) ' 

d|q| , 
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D n ~ D22 



dp 



2 q + P 



f x (p)q p e° x » K f v (p)(2q + p) a 



x<j (hah/3 + hahp) 



2ei K ei p /^ 1 1 ) (q) - 2e 2K e 2p , / 9 2 1 2 ) (q) + (e 2K e lfl ~ e lK e 2fi )[p[^ (q) + /0^(q)] 
2e 1/3 ei CT pf 1 ) (q + p) + 2i 2 pha P22 (<1 + P) 

+(ei/3e 2(T + e 2 /5?i CT )[pf 2 ) (q + P) + /4i (q + P)] 



Off 



2 ^ (0) (|qi: 



D 



12 



Dai = ~(^ff) 5 



dp 



2|q + P 



TA(p)g P e CTApK T ; ,(p)(2 9 + p) Q 



x< (ei CT ei/3 + £20-62/3) 



(ei/teip + e2 K e2 M )[pi 2 ) (q) + Pzifa)] + ^l^Pui^) + ^2^1^ (q) 



d|q| 



Ji), 



■(e2«ei/i + £ik£2m) 



2ei/3£i CT pl 1 1 ) (q + P) + 2? 2/3 e 2ff p 22 ) (q + P) 
+ (ei/3£2 CT + £1^2/3) [P12 (q + p) + P21 (l + P)] 

Di2-D 2 i = -(2zg) 2 / ^ f A (p) gp e CTA ^(p)(2 g +p) a 

j * q 1 pi 

x< {hah/3 + hah/3){eiK(-ifi + e 2K e 2M )[pi 2 ) (q) - P21 (q)] 



Off 



, 2 ^ (0) (|q|) 



d|q| 



-(e2 K ei M - ei K e2 Al )(ei^e2<r - eio-e2^)[Pi2 (l + P) ~ P21 (Q + P)l | • 



Given the above expressions it is easy to see that no mixing occurs between Ay and the set {Aj, Aq, Ajj} as a result 
of the damping term Eq. (|43[) . Therefore, in our approximation, no circular polarization is generated by first order 
processes up to 0(g 2 ). 



V. THE SECOND-ORDER INTERACTION TERM 

We now move on to calculate the contribution to the evolution of the photon density matrix from scattering 
processes which are second order in the interaction Hamiltonian density operator Eq. ([34]) . Details of the calculation 
of the second-order interaction Hamiltonian and the corresponding refractive term are presented in Appendix [Bj The 
second-order interaction Hamiltonian is given by 



2 id int (i) = -^(2gf 



((Pi +P2 + hU2 P 2 + J 2 ) x ) e -(b'O o -0O o -(P») u -(W o )t tf 3 (pi _ lx _ p2 _ h) 



(51) 

(P2)e s /3 (pi )at (P2)a r (pi ) 
(pa + h) 2 + ie 



, It , 1 \ to j \ \ i((p 1 )"+h ~(P2) +h )tx3t 1 ,i , e *sf3(Pi) e r K (P2)ai(pi)a r (p 2 ) 

+ {(pi +P2- h) a (2p2 - h)x)e l((Pl> +tl {P2) + ' 2 }t 5 A pi +Ii - p 2 + 1 2 )— 7 , , 9 , . 

(pa - '2) + «e 



The refractive term of Eq. ([27)1 due to this interaction Hamiltonian is 

= i(27r) 3 S 3 (0)(2g) 2 (5 urPsv (q) - Wm-(q)K(9K(<?)V(g), 



2 fli„t(0),P u „(q) 
where we have defined %j(q) via 



(2 7 r) 3 5 3 (0)7-,( (? ) 



M2 (2tt) 3 ( 5 3 (1 1 + la)^, 1 2 )T < (Z? ) \ 2 )[q ■ (h + k)} 

1 1 

(Z2) 2 + 2l 2 ■ q + ie ~ (l 2 ) 2 -2l 2 -q + i€ 



(52) 
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Note that the interaction Hamiltonian Eq. (|5T|) and the refractive term Eq. (|52f were computed using the photon 
propagator in the Feynman gauge Eq. (|B-3[) , but a demonstration of gauge invariance of these results is detailed in 
Appendix [C] For an unpolarized photon ensemble, it is easy to see that the right side of Eq. ((52|) vanishes when 
evaluated in terms of p^ 1 defined in Eq. above. The contribution to the perturbed density matrix, p^\ to second 
order in the interaction coupling g is 

|p$(q) = ^{5urpi 1 J(ci)-S vs p^(cO)e^ r ( q )T^ q ). (53) 

Here we have used Eq. (f27|) and have ignored the damping term, which is of 0(g 4 ). Expressing Eq. (f53|) explicitly in 
terms of photon density matrix components gives the evolution equations 

|p<i>(q) = ^(-^^(q) + ^^(€0)^(5) (54) 

|p&(q) = ^(^^(qj-e^^q))^^) (55) 

|pW(q) = ^£(tfe? - «/$(q) - e^ft (q) - p£(q)]) T^fo) (56) 

|p«(q) = ( - ~ ^e^fa) + efe^ftq) " P^Cq)]) V(<?)- (57) 

All polarization vectors e(f in the above expression depend on q, the same photon momentum as in the argument 
of the photon density matrix. Using the Stokes brightness perturbations in Eqs. (|30[) - (|3"3"|) , we can now express the 
evolution of the polarization of the photon ensemble due to processes mediated by 2 Hint (t) as 

|Az(q) = 0, (58) 

|A Q (q) = -^(C^Ay^+z^A^q)), (59) 



^(q) = -^(x(q)A v (q)-^(q)A (q)), (00) 



|A v (q) - ^(x(q)A c/ (q) + C(q)A Q (q)), (61) 
where we have defined the contractions 

C(q) = -2V(e^ + £ ^), (62) 

X (q) = 2T^(e^-e^ 2 ) (63) 

V(q) = 2T p „(e^-e^) (64) 

and the quantity 7^„ is defined via the integral expression Eq. (|52p . 

The circular polarization brightness Ay is sourced by terms which are proportional to either Aq or A^: a linearly 
polarized photon ensemble in the presence of the interaction Eq. (|34[) will acquire circular polarization due to processes 
which are of second order in the interaction. As long as the interaction acts, both Stokes brightnesses Aq and Ajj are 
rotated with Ay. Note that the above equations do not depend on the time component of T^(x) since the polarization 
vectors are purely spatial. 



VI. DISCUSSION 



The evolution equations ([55]) to (fuTj) , along with Eqs. (|4"6"|) to (j4"9"|) , are the central result of this paper. Other 
source terms associated with the usual Compton scattering effects will also appear on the right sides. While the 
polarization brightnesses will be zero prior to recombination, during recombination Aq and Ajj become non-zero, 
with an amplitude a factor of 20 smaller than the intensity fluctuations A/. It is easy to see by inspection of the 
evolution equations that at that point, Aq and Ay will rotate into each other with a characteristic angular frequency 
ujqv — g 2 (/k°, Au and Ay will rotate into each other with a characteristic frequency uyy — g 2 x/k°, and Aq and 
Ajj will rotate into each other with a characteristic frequency ujqjj = ga, along with an exponential decay or growth 
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of Aq and Ajj associated with the first order damping effects. All of these source terms are active whenever the 
interaction Eq. (|34[) is nonzero, in contrast to the conventional Compton scattering terms, which are only significant 
when the photons propagate through ionized regions of the universe. 

The rotation between Aq and Ajj can be constrained from current measurements in a straightforward way. Linear 
polarization on the sky is conveniently exp ressed in a different basis, corresponding to the "gradient" and "curl" 
pieces of the polarization tensor field [3y, [37[, also known as the E/B decomposition [33|. This decomposition is useful 
because scalar density perturbations in the universe, which evolve into the structures we see today via gravitational 
instability, generate only E-mode polarization. Subsequent rotation of the polarization plane as the wave propagates 
rotates E-mode into B-mode. Current limits on the amplitude of B-mode polarization (see Refs. [H, H(| El for some 
recent linear polarization measurements) can be translated into limits on the total rotation of linear polarization 
between the time of last scattering and today; see Refs. [2(], EH for corresponding limits on magnetic fields due to 
Faraday rotation. Precise limits on the interaction studied here from the first-order rotation effect Eqs. (|4"T)) and (j4"8|) 
can be obtained similarly, and will be computed elsewhere. But we know that the total amount of rotation must be 
small at frequencies between 50 GHz and 150 GHz where good measurements of the primordial linear polarization 
have been made. Given this observational constraint on linear polarization rotation, can some realistic cosmological 
model generate detectable circular polarization via Eqs. (|59 p - (j6"Tj) ? First, note that u>qu has dimensions of [gT], and 
that for this rotation to be below current limits, 

uj QU ~ gT < H ; (65) 

otherwise, as the microwave background photons propagate a Hubble distance from last scattering until today, we 
would have substantial rotation of E-mode into B-mode polarization. We are aware of no other laboratory or theoretical 
constraints on this class of interactions. 

Now note that any first order damping effects will contribute an extra factor of length compared to a (q), arising 
from an additional time integral over the field T^(x); if the field is active for all times, this leads to roughly a factor 
of Hq . Then the time scale for exponential growth or decay of linear polarization is, by dimensional estimate, 
u)qu(u)qu I Ho), which is small compared to loqu- we can always neglect the exponential growth or decay of linear 
polarization compared to its rotation. 

For generation of V polarization, loqv and couv both have dimensions [g 2 T 2 /k ], differing only by a geometric factor 
related to the propagation directions of the photons and the direction of the field T^(x). So a dimensional estimate 
for both is ujqv — ^uv — ^Qui^QU )• A typical microwave background photon today will have a frequency of 
k° ~ 100 GHz or 10 11 s _1 , while the characteristic size of ujqu at the observational limit is Hq, a huge mismatch in 
scales. So in the case considered here, the generation of circular polarization is always vastly subdominant to the 
rotation of the linear polarization, and can be neglected. 

The calculation presented here demonstrates generally that, given additional interactions beyond Compton scat- 
tering, circular polarization is not necessarily zero, and elaborates the framework for calculating it for a given micro- 
physical interaction. Other interactions may well induce circular polarization without optical activity from the linear 
interaction term, and for these cases circular polarization could be the most constraining probe. We have also only 
considered a constant field T p for simplicity; calculations for a non-constant field are messier but straightforward, 
involving convolutions over the field and photon distributions. Spatial or temporal variations in the field could change 
the relative importance of the optical activity and circular polarization generation effects. In particular, the torsion 
field necessarily couples to fermions via the interaction Ctf — gi1~ fJ „ipJ >1 '4> as well as a torsion-induced four-fermion 
interaction Cff = 52V>7 5 7' j V'V'757mV' OS, EH], where g\ and g\ are renormalized couplings. It would be interesting 
to study the coupled system of torsion and fermions subject to our formalism for evaluating V, since backreaction 
effects from the torsion-fermion interaction could enhance the amplitude and modify the scale of V; we leave this to 
future work. 

We encourage experimenters to make measurements testing the standard lore that circular polarization of the cosmic 
microwave background radiation should be identically zero, and theorists to consider the effects of any non-standard 
photon couplings on microwave background polarization as photons propagate over cosmological distances. 
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APPENDIX A: FIRST-ORDER CALCULATION 



In this appendix we determine the contribution to the time evolution of the photon density matrix due to processes 
which are first order in the interaction Hamiltonian density Eq. ()34[) . We first detail the calculation of the refractive 
term and then the damping term of Eq. (|27|l . The first-order interaction Hamiltonian is simply given by 



(A-l) 



In the approximation employed here, all operators in the collision terms of Eq. (|27|) are from the free theory. From 
this point on, all operators represent free theory operators and we drop the '0' superscript used in Sec lIIII Ignoring 
any processes in which two physical photons are either annihilated or created leads to the following expression for the 
interaction Hamiltonian: 

1 H int (t) = 2g J d 3 xdpdkdk' e"" " f v {p){-ik' a ) 

x (at(fc)a r (fc')e: Al (fc)e^(fc')e l ( f,+fc - fc ')- - a\. (tf)a.(k)ZpW)e a ^)e^ k '- k > x ) , (A-2) 

where we have used the expression for given in Eq. (|14p and we have used the shorthand notation 



dp : 



d A p 
(2^ 



dk 



d 3 k 



(27r) 3 2fc° 



Performing the spatial integral and relabeling dummy momentum variables and spacetime indices in the second term 
of (|A-2|) gives 



l Htat(t) = -2ig J dpdkdk\2^f8 {i \k + V -k')e t{ko+Pa - k '« )t e^ aP f u {p)a) s {^ (A-3) 
Now perform the J dk' integral, after which we have 



l H int (t) = -2ig I dpdk ^ j 2|k + p| ot(fc)a r (fc)e; M (fc)(^(|k| + |k + p|) + ^(2fe i +p i ))er/j(A), 

— n t 



2ig / dpdk- 



2|k + p| 



-dl(k)a r (k)e* (k)(k + k) a e r/3 (k) . 



(A-4) 



where (k) a — (|k + p| , k + p). Equation (j A-4j) is our first-order interaction Hamiltonian. Now the commutator 
necessary for the refractive term of Eq. (|27p is given by 



l H int {t),V uv {q) 



2ig / dpdk 



f(p) e i(ka+Po-\k+p\)t 



(k)(e^(k + k) a )e rfj (~k) 



2|k + p| a " 
:(2Tr) 3 2q° (<5 ur 5 3 (q - k - p)at(k)a v (q) - S vs 8 3 (q - k)al(q)a r (k) 



(A-5) 



where we have used the canonical commutation relations between the free creation and annihilation operators Eq. (|15p . 
Taking the expectation value of (|A-5I and using the relationship between the number operator and density matrix 
given by Eq. (|17|) . we arrive at the following expression: 



Perform the J dk integral gives 



29 J dpdk 2|k + p| e^(k)e^(k + k) a e r0 (k) 

x(2n) 6 (2q°) 2 5 3 (q - k - p)S 3 (q - k) (5» r p st) (q) - <W«r(q)) . 



(A-6) 



^MV^ia) ) = 2. 9 / dp "^f. , , ■ , e* (q)e^(q + q) a e r0 (q) 



4|q||q + p| ^ 

x(2n) 3 (2q°) 2 S 3 (~p) (6 urPsv {q) - S vs p ur (q)) , 



(A-7) 
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where as before we have defined (q) a = (|q + p| , q + p) so that the above can be expressed as 

T , j / (p)e% 0+Po- ' q+p D t 



^(t), £«,(<!)]) = 2ge* stl (q)(2q°)(6 ur p sv (q)-6 vs p ur (q)) j dp (2^) 3 ,5 3 (-p) tJL ^ ■_ _ _ ' ' ' (A-8) 
x (e^ 0/? (|q| + |q + p|) + e^(2 ?j + Pj ))e r/3 (q + p) . 

Define the quantity 

(2tt) 3 £ 3 (0)A^(q) = (2^) 3 5 3 (0)^ [M/3] (g), 

= | ^(2^ 3 (-p)^^(^(|q| + |q + p|) + ^'' 3 (2g J +^)), (A-9) 



where we have anticipated the fact that the delta-distribution factor will be present once an appropriate form for the 
external field T^x) is chosen. The refractive term of Eq. ([27|) due to first order processes can then be expressed as 



^intCO), A„(q)J ) = 4g g °(2^) 3 ( 5 3 (0)e: AI (q) (6 ur p sv (q) - S vs p ur (q)) (q)e rP {q) . (A-10) 
Next we consider the damping term of Eq. (|26p , the integrand of which involves the double commutator 

J' l/ rp) e »((*i)a+Po-|kl+p|)t f x ^ e i((k 2 )o+la-\k 2 +l\)(t~t') 



1 H- mt (t-t'), 1 H- mt {t),V uv {q) 



(2igY / dldpdk 2 dk 1 



2|ki + p| 



2 |k 2 + 1| 



(A-ll) 



xe s * /1 (fc 1 )e^(fc 1 + fci) Q e r/3 (fc 1 )e; CT (fc 2 )e CTA " K (fc 2 + k 2 ) p t rnK {k 2 ) 

x(2 7 r) 6 2< z ° (s ar 6 3 (q - k x - p)(2(fc 1 ) ( 5 ms ,5 3 (k 2 + 1 - ki)fi+ (fc 2 )a„( g ) 

-2g°<5 m ,£ 3 (k 2 -q)aj(fci)a m (fc 2 )) 
-^ s <5 3 (q - ki)(2<j°(5„ m( 5 3 (k 2 + 1 - q)al(k 2 )a r (k) 

-2(k 2 )°6 nr 6 3 (k 2 - k x - p)fi£(9)Mfe)) 

where we have used the canonical commutation relations between the free creation and annihilation operators, Eq. (|15[) . 
and have defined (fci) Q = (|ki + p| , ki + p) and (k 2 ) a = (|k 2 + 1| , k 2 + 1). Now the J dki and the J dk 2 integrals can 
be performed, giving 



1 H il *(t - t'), 'HtatW.^^q) - (2 lff ) 2 2< Z ° / dldpf x {l)e° x ^f u {p)^ afi 



(A-12) 



4|q-p| |q| 
e t(|q-p-i|+io-|q-p|)(t-t') 



2|q-p| 



4 q-p-1 |q-p| 



e *na{q-p-i){q-p-i + q"-p)pem K (q-p)8msat(q-p-i)a v (q) 



-2q a 



i(«o+Jo-|q+l|)(t-t') 



4|q||q + l| 

e <((9)o+PO-|q+p|)t e ^( 9f )(g + g +p ) a e r/3 (g +p ),5, 



(9+/) 



2g° 



Cr(<7-0(<7-i + 9) P <w(g)6 m „a) 1 (9-/K(<7+p) 



4|q + p| |q| 

e t(|q-l|+J -|q|)(t-t') 

4|q-l||q| 

e t(|q+p|+io-|q+P+l|)(t-t') 

(q +p+ l)6 

nr 



4|q + p||q + P + l| 
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After taking the expectation value, this becomes 

'^^(t-t'), \ 1 H i nt(t),T> uv (q)]]) = (2igf I dldpf x (l)e^f u (p)e^ a f } (2n) 3 S 3 (l + p)e^ + ^ t (A-13) 



2|q + l| 

(q) - e n CT (9)(? + q+l) p e mK (q+l)SnvPsm(<l + 1) 
e-igq+pl+io-lqlK^fr + g +p ) aC ^(g +p )^ a 

2|q + p| 

e na(<l+p)(q+P + q)p^m K (q)S m uPnr((l + P) - e£ a (q+p) (q+p + q) pe mK (q)Snr Pum(q) 



Interchanging Z and p in the first set of terms gives 



^(t-O. [^intW.^^q)]]} = (2igf J dldp{2ir) 3 6 3 (l + p)t 



,i(po+lo)t 



e- l(|q|+P0 - |q+p|)t VA^(rte gAp ^^(/)£^" /3 e: At (g + p)(g+ P + g)a£rMg)^r 

2|q + p| 

<<t(9)(9 + 9+p)pW((?+p)<*m S Pni;(q) ~ fa) (<? + g+p) p em«(<Z+jp)<W S m(q + P) 

e-'(lq+Pl+'°-|ql) t, f A (i)e g ^f y (p)e^e; M (g)(g + g +p ) a e^(g +p )^ a 

2|q + p| 

£*na(q+p)(q+P + q) P £ mK (q)$ m uPnr(<l + p) - e* CT (<? +p ) (q+ p + q) p e mK (q)5 nr p um (q) 

After relabeling some spacetime and polarization indices in first set of terms, the above becomes 

^H^t-t'), [ 1 H int (t),V U v(ci)\]) = (2igf J dldp{2^) 3 5 3 (\ + p)e^+ 1 ^ 
fUp)e^ K f x (l)e^e; ia (q +p )(q +p + q) p e mK (q)e* all (q)(q + q +p ) a e r0 (q +p ) 



x | e -*(|q|+Po-|q+p|)t'^ 



.p-idq+Pl+io-lql)*'^ 



2|q + p| 

OrnPsvyQ) ~ ®svPnr (q + p) 

(q + p) 



} 



Now integrate over dt' (t — > oo) and define Ag = |q + p| — |q| to arrive at 



/OO 
-oo 

(2zg) 2 | dZ dp jqj 1 ^ P) r y (p)c^f A (0e^e ng + g)„(g + g +p ) 



x<U(p - Ag) 



^rn^um^mfc^s/iPsv (l) ^um^sv^mK^sfxPrir (O. H~ P) 



+<5(Z + Ag) 



(q + p) 



}• 



(A-14) 



(A-15) 



(A-16) 



where for convenience we have defined D uv above, as well as the abbreviations e rfJ _ = e r/J (q) and i rfl — e rll (q + p). 
Note that since p is expressed in a linear polarization basis all polarization vectors above have been assumed to be 
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real. Now in the J dl integral above, the relevant factor can be simplified as 

p [ dlS(l + p)f x (l)(q + q +p ) p e° X ^f„(p)(q + q +p ) a e^ 



X 



= e^w^lf^o, " p )[l q l + h + PW MK - T a (k, ~P)(H + Pi)*?* * 

T>o,p)[|q| + |q + p| -p ]e^ - 2T ( Po , p)q^ 
= e„ CT e mK e^e r/3 |f i (Z ,p)[(|q| + |q + p|) + lole™ " + 2T (l , p)q l e al ° K 

x jr> ,p)[|q| + |q + p| -p ]e^ - 2f (p ,p) qj e^ 

where in getting to the final expression we have used pwr v \ (p) = 0, the fact that e np (k) can be fixed as purely spatial 
for any momentum k, and the assumed behavior of T p (x) under a parity transformation, Eq. Q36p . Now perform the 
appropriate energy integrals made trivial by the presence of delta-distributions in Eq. (|A-16|) . The relevant factor 
from the second term of Eq. (|A-16|1 becomes 

e n(T e mK e sfi e r p J dl S{l + Ag) jf^o, p)[|q| + |q + p| + h]e ai ° K + 2T (Z , p)^° K } 

x|f i ( Fo ,p)[|q| + |q + p| -Po]^ 00 -2f (p a ,p)q^ Q ^ , 

= e na e mK e^e rP S.f i (-Aq,p)[\q\ + |q + p| - Aq]e° i0K + 2T (-A (? , p)^e CTfflre J 

x|f,(p ,p)[|q| + |q + P | -p*V m - 2f (p ,p) qj e^ 
= e„ CT e mK e s/J e r J - 2T l {Aq,p) |q| e at ° K + 2f (Aq 7 p)q l e' TtOK 

x|f,(p ,p)[|q| + |q + P | ~Po]e^ - 2T (p , p)<fce w ' 0/3 } , 
= e na e mK e sl ,e r p( - 2T x {Aq, p)q P e aXpK ~^ jf^po, p)[|q| + |q + p| - p ]e^ ~ 2f (p ,p) qj e 



(A-17) 



where in arriving at the final expression above we have used the assumed behavior of T^x) under a time reversal 
transformation Eq. (|37p ■ Now perform the po integral present in the first term of Eq. (|A-16[) and arrive at 



inae mK e Sil € r p J d Po S{p - Ag) |fi(Z 0) p) [|q| + |q + p| + IqK Mk + 2f (l ,p)q i e' TU 

fj(po,p)[\q\ + |q + P | - P o]e MOP - 2f (p ,p)q 3 e^ 
e w e mK e sll e r pi Tj(Z 0) p)[|q| + |q+ p| + lo}e aiOK + 2T (l , p)^ "}/ 2^?, p)q a e^ 



Next change the integration variable Iq — > — po remaining in the first term of Eq. (|A-16p and the factor above becomes 

e nrT e mK e s „i r pS. - f 4 (p , p)[|q| + |q + p| - Po]e CTl ° K + 2T (p Q , p)q i e rTi ° K \ Uf u (Aq, p)q a e^\ , (A-18) 

where we have again used Eq. (|3"T|) . Following some relabeling of spacetime and polarization indices, the factor (|A-18[) 
becomes identical to the factor (|A-17|) . Thus Eq. (|A-16|) can be expressed as 

D uv = -{2ig) 2 e mK e sp J f x (p)q p e aXpK f v (p)(2q + p) a e^l n<y l r p (A-19) 



x S nr 5 us p mv (q) + 5 rn 5 vs p um (q) - S us 5 mv p rn (ci + p) - 5 um S sv p nr (q + p) 
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where we have defined the 4— vector p = (Aq, p) and recall that Aq = |q + p| — |q|. The linear combinations necessary 
to describe the evolution of the independent polarization degrees of freedom are 



D n ±D 22 = -{2igf 



dp 



2|q + p| 



T x (p)q p e aX P K T v { P ){2q + p) t 



(A-20) 



x<^ (hvhp + £2(762/3) 



2ei re ei jU/ on(q) ± 2e 2K e 2AI ,c>22(q) 
+ (e 2K ei AI ± ei tt £2 M )[ ( oi2(q) +P2i(q)] 



£>12 ± D 2 1 



dp 



2 q + P 



2ei/3£i CT pii(q + p) + 2e 2/5 e 2(T/ 9 22 (q + p) 

+ {.h(tha + e2^ei<r)[Pi2(q + p) + P2i(q + p)] 
TA(p) gp e CTA " K f,(p)(2q+p) Q 

eiK£2 P [pii(q) ± pn(q)] + £2 K ei M [p22(q) ± P22(q)] 
+(ei K ei A1 + £2 K £2 M )[pi2(q) ± P2i(q 



(A-21) 



If we define the quantity 

(27r) 3 ^(0)£ng) 
and the integral operator 

(27T) 3 S^(0)B 2 [q;p. 



ei/3eio-[pn(q + P) ± Pn(q + P)] + e2/3e2<r[/0 2 2(q + p) ± P22(q + P) 
+ (€ ll3 e 2a ± hahp) [pi2(q + p) ± P2i(q + p)] 



dp 



2|q + P 



f x (p)q p e' 7X P K f l/ (p)(2q + p) a e^ e na 5 nr e rl3 



uv\ — €rriK,€sfi 



dp 



_ -T x {p)q p e^f„{p)(2q+p) a e^i na i r p 
^ iq + PI 

X (5u S ^ m «yO™(q + P) + <5 U m^ S i,Pnr(q + P)) , 



(A-22) 



(A-23) 



Eq. (|A-19j) can be expressed in the general form 

D U v = -(2«£j) 2 (2tt) 3 6 (3) (0K e mK (e UAt /w(q) + e VfJ ,p um (q)) Bf K (q) - B 2 [q; p U v(q)] 



(A-24) 



APPENDIX B: SECOND-ORDER CALCULATION 

In this Appendix we explicitly calculate the relevant quantities describing the evolution of the photon density matrix 
due to processes which are second order in the interaction Hamiltonian Eq. (|34p . The second-order scattering matrix 
operator is 



5(2) = _i 



OO pOO 



dt 



OO J — OO 



dt'T{H^(t)H^(t')}, 



(B-l) 



where H; nt (t) is the first-order interaction Hamiltonian, Eq. (|A-1|1 . We will denote the interaction Hamiltonian 
operator which has a non-zero overlap with a single photon lying in both the initial and final scattering states as 
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2 H mt (t). Applying Wick's theorem to simplify the time ordered product in Eq. (|B-1[) gives (ignoring vacuum terms) 



^ 2 \t) = ~(2gf f d^ye^e^T^T^y) : fA^x)d a Ap(x)A p (y)d x A K (y) 



+A fl {x)d a A (x)A p (y)d x A K (y) + A M (x)d a A p {x)A p {y)dxA K {y) + A^(x)d a A (x)A p (y)d x A K (y) j : , (B-2) 

where all partial derivatives are understood as acting solely on the function immediately to the right; the variable 
being differentiated is in the argument of this function. We denote the contraction of two operators A and B by 

n 
AB . 

To simplify Eq. (jB-2|) . we use the free-theory photon propagator in the Feynman gauge, 



D^(x — y) = 



d 4 k -ig^e- ik < x -^ 



(2tt)< 



k 2 



(B-3) 



(See Appendix JC]) below for a demonstration of gauge invariance, where we explicitly consider a different gauge-fixed 
photon propagator). In order to deal with the derivative couplings, we interpret the time ordering as T* ordering; 
specifically, we require that derivative couplings act outside of the time ordering operation [35[. For convenience, 
define the operators 



A+(x,p) = a s (p)e sf ,(p)e ip ' x , 
A-(x,p) = at(p)e* (p)e*-*. 



(B-4) 
(B-5) 



Equation (|B-2j) then becomes 



iH< r '(t) 



d 3 xd 4 y 



d A k d 3 pi d 3 p 2 



T v {x)T a (y) 



-ic 



-ik-(x-y) 



IC 



(2tt) 4 (27r) 3 2p° (27r) 3 2p0 

g l _ ip {-ipi a ){A+{x,pi) - Ap(x,p 1 ))(-ip 2 \)(A+(y,p 2 ) - A~(y,p 2 )) 

+ (ikx)g^ K (-ip la )(A+(x, Pl ) - Aj 3 (x, Pl )) (A+ (y,p 2 ) + A~ (y,p 2 )) (B-6) 
+ {-ik a )gpp(A+(x,p 1 ) + A-(x,p 1 ))(-ip 2 \)(A+(y,p 2 ) - A~(y,p 2 )) 

+ (-ik a )(ik x )90 K (A+(x, Pl ) + A-(x, Pl ))(A+{y,p 2 ) + A-(y,p 2 )) 

Picking out the non- vanishing overlap of H ^ (t) on single-photon initial and final scattering states and calling this 
2 H int (t) gives 



l H int (t) = ~2 (25)2 / I dkd Vl d V2 e^eP^T v {x)T a {y) 



-te 



-ik-{x-y) 



k 2 



g IJ , P {-ipia)(-ip2x)( - A K (y 7 p 2 )A+(x, Pl ) - A fj (x, Pl )A+(y,p 2 )) 

+ 9»k (ik x )(-i Pla )(A- (y, p 2 )A+(x, Pl )- A~ {x lPl )A+ {y, p 2 )) (B-7) 
+ 9p P {-ik a ){-ip 2 x)( - A~(y,p 2 )A+(x,p 1 ) + A~(x,p{)A+(y,p 2 )) 

+ gp K (-ik a )(ik\)(A- (y, p 2 )A+ (x, pi) + A~ (x, pi)A+ (y, p 2 )) 

where we have again used the shorthand notation for the integral measures defined in Appendix (jA|) . After some 
relabeling of spacetime indices, Eq. (|B-7j) becomes 

x {-ig» P ) ( ° u2 , Tr ) (pia(P2 - k) x - k a (p 2 - k) x )A~(y,p 2 )A^(x,p 1 ) (B-8) 



! H int (t) = --(2gY J d a xd*y J dkd Pl dp 2 e>~ < 

-ik-(x-y) 



k 2 + ie 

e -ik-(y-x) 

k 2 + ie 



(pi a (p2 ~ k)x - k a (p 2 - k)\)Ap (x,pi)A+(y,p 2 ) 
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using the property D fiu (x — y) = D fJiU (y — x). Now plug in the Eqs. (|B-4p and (IB-5|) . express T^(x) in terms of its 
Fourier transform, and perform the J <ix and J dy integrals to get 



2 H int (t) = - l ^L{2gf I dkd Pl dp 2 



e^tP aXK f v {h)f a {l 2 ) 



k 2 



(-ig t i P )(pia{p2 - k)\ - k a (p 2 - k) x ) 
i(( P1 ) u -(iO u +* u )* tf 3 (pi _ h + k)s i {p2 + h + k)e*Me rfi {p x )al{j> 2 )a r { Pl ) (B-9) 

Next perform the k integral: 

2 H int (t) = -!^( 2ff ) 2 y' ( ip 1 dp 2 ^ 1 dZ 2 e^e^ A ' t f i ,(/ 1 )f (7 (; 2 )(- i5/ip ) (B-10) 

((pi +P2 + h) a (2 P2 + fej^e-cw-^-w-o-^oP! u P2 i 2) ^^);^;)g^)^^) 

(p 2 + t 2 ) +?e 

+ ((pi + P2 - ^ 2 ) Q (2 P2 - fa) A )e^ W,.,.^ + h _ p2 + l2) ^^K.^)at( Pl )^fe) 

(p 2 -/ 2 ) 2 +ze 

We are now in the position to compute the commutator of 2 H- m t(t) and 2? uv (q) necessary for the refractive term of 
Eq. (p7) . This is given by 



! ff int (t),I?™(q) 



z(2tt) 



(2<?) 2 / dpidp 2 dl 1 dl 2 e> iVaf3 e>" TXK T l ,(h)T a (l 2 )(-ig flp )e 



(B-ll) 



((Pi +P 2 + y Q (2p 2 + / 2 ) A )e-^)°-^)°)^ 3 (Pi - li - P2 - h ff^§^ 

(P2 + h) + ze 

x (27r) 3 2g° (S 3 (q - Pi)8 ur d\{ P2 )a v {q) - 5 3 {q - p 2 )S vs al(q)a r ( Pl )) 

+ ((px +P 2 - ; 2 ).(2p 2 - Z 2 )A)e i «^)°-^)°)M 3 (p 1 + U - p 2 + 1 2 )^^M 

(P2 - <2j + «e 

x(2tt) 3 2<j ((5 3 (q - p 2 )5 ur aj(pi)a4<?) - <5 3 (q - Pi)§v S al(q)a r (p 2 )) 
Take the expectation value of Eq. (|B-11[) above to get 



l H int (t),V uv (q) 



z(2tt) 3 



(2g) 2 / dp^dhdhe^e^fvih^ih^-ig^e 



i{h°+l2°)t 



(B-12) 



x(2vr) 6 (2( ? ) 2 ( 5 3 (q - Pi)<5 3 (p 2 - q)<5 3 (pi + I x - p 2 + 1 2 ) (5 ur p sv (q) ~ S vs p ur (q)) 

x \ P1 + P2 + hU2 P 2 + hhe^r-toryt &M^M 

{P2 + hY + «e 



Hp, +P2- hU2 P 2 hhe^ - ^'^^) 

VP2 - hY + if 



Now perform the pi and p 2 integrals, giving 

Z(27T) 3 



■H int (t),V uv (q) 



urPsv \Q) OvsPur (q)) 



(2g) 2 / d/ 1 d/ 2 <5 3 (l 1 + l 2 ) e ^ e " CTA ' i T' l/ (Z 1 )f' CT (Z 2 )(- J5/ip )e l ( /l +< 2 >* (B-13) 



(2g + l 2 ) a (2q + b)A ^#^ + (2. - fa)^ - * 2 )>> (<ZKk(<Z) 



( g + ? 2 ) 2 + ie v J zy " v * z/A (<z-Z 2 ) 2 + *e_ 
Focus now on the following factors present in Eq. (|B-13|) : 

e^e^ XK f v {h)f a {l 2 )e* SK {q)erp{q){2q ± l 2 ) a (2q ± ^ 2 ) A = 2e^e/ AK f i ,(; 1 )T' (T (; 2 )e: K (g) e ^(g)(2g ± / 2 ) q(?a . 
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where we have used InT^l) = in arriving at the expression on the right. Now use the epsilon identity i^ va ^ t^aXn 
— 6(5^(5? 5$ to perform the contractions above and arrive at 



-2 ( [Ti • T 2 ] [{2q ± h) ■ q] [e r ■ e s ] - [T 2 ■ T x ] [e r ■ q] [(2q ± l 2 ) ■ e s ] + [T 2 • e r ] [(2q ± l 2 ) ■ e s ] [T x ■ q] ( 13-1.4 i 

-[f 2 • e r }[(2q ± l 2 ) ■ q][e, • f x ] + [(2<j ± J 2 ) ■ f 2 ][q • e r ][7\ • e s ] - [f 2 • (2<? ± I 2 )][fi • g][e P • e,] 



To simplify this, use e r ■ e 8 — ~S rs , which when contracted with (S ur p sv (q) — S vs p ur (q)) in Eq. (IB-13j) vanishes, as 
well as e r (q) ■ q = and q ■ q = 0. The remaining terms from the expression (|B-14|) then become 



T-2[f (ll 1 2 ) • e r ]e sll q v (^2)^°, ~h) - {hY?^ -l a ) 

= ±2[f (& la) ■ £ ,]e^o f (h) l f °(f°, 1 2 ) + (ia) ^^, 1 2 ) ) = ±2[T(i°, 1 2 ) • e r ]e si f l 2 )g (Z° + 1°) , 



where we have used Eq. as well as repeated use of Z[aT,t](0 = 0. Inserting this expression back into Eq. (|B- 131) . 
we have 

" H int (t),V U v(q)] ) = -(2n) 3 (2g) 2 (d ur pM - 8 vsPur (q)) J dhdl 2 6 3 (h + \ 2 ) e l ^°+ l ^ t (B-15) 

1 1 



x[ £r .fG 2 °,l 2 )][f(/0,l 2 ).e s ][(7-(/i+; 2 )] 
We now define a quantity 7L- by 



(l 2 ) 2 + 2l 2 -q + ie (l 2 ) 2 -2l 2 -q + ie 



dhdh (27r) d ^(li + h)Tj(ll h)Ti(l», l 2 )[q ■ (h + h)} 
1 1 



(B-16) 



_{h) 2 + 2l 2 -q + ie {l 2 ) 2 ~2l 2 -q + ie_ 
Equation (|B-15|) can be expressed conveniently in terms of this quantity: 

r2 £int(0),A„(q)l) = i{2ir) 3 5 3 (0)(2g) 2 (S urPsv (<l) - W«r(q)K (<?)<(<?) : 



(B-17) 



where we have exploited the fact that e l s (q)el(q)Tij(q) — e^{q)e^{q)T )JiV {q) since e(q) is purely spatial. Explicitly, the 
components of i 



l H int (t),V uv (q) J are given by 

^int(*),Ai(q)]) = *(2^ 3 <5 3 (0)(2. 9 )% l/ ( (Z )(-65 1 (g) e ^(g)p 12 (q)+e^( 9 )e5'(g)p 21 (q)), (B-18) 
J ffi„t(i),Pi 2 (q) ) = *(2^) 3 ,5 3 (0)(2. 9 ) 2 V(g)([ e 5 1 (g) e r(g)- e ^g) e ^g)]p 12 ( q ) (B-19) 



, l ( ? )e^(g)[ Pll (q)-p22(q)]), 



(B-20) 



2 ^intW^2i(q)J) = i(2n) 3 5 3 (0)(2g)% u (q)(- \4{qK{q) - e v 2 {q)e»{q)]p 2l {q) 

+e?(qy 2 (q)[p 11 (q) - p 22 (q)}) , 
2 ffint(t),P 22 (q)l ) - »(2^ 3 J 3 (0)(2. 9 ) 2 V( g )(ei I (q)e^(q) Pl2 ( q ) - e^q)e^q)p 21 (q)) . (B-21) 



APPENDIX C: GAUGE INVARIANCE 

Here we verify that the calculation in Appendix B is gauge invariant by explicitly using a different gauge-fixed 
photon propagator, namely 



d 4 k —ie~ ik '( x ~ y *) 
D, v {x-y)= I (2?r)4 k2+lt [ft- -(1-0 



(C-l) 



19 



Isolating the contribution to 

2 ff int (i) c ,P„ t ,(q) 
x (S U rPsv(q) - S vs p ur (q)) 



2 Hi nt (t),f> uv (q) ) due to the term linear in (1 — £) gives 



(q + h)^(q + h) P , , , , n / _, i t \ e* SK {q)e s p(q) 
^ + ^ 2 (q + q + h) a {q + q + h)x 

{q-k)^q-l2)p 
{q-hY 



(q + q- h)a(q + q~ h)\ 



(q + h) 2 + ie 
t* sl 3(q)tr K (q) 



(C-2) 



(q-h) 2 + i( 



where we have performed all the integrals similar to those in arriving at Eq. (|B-13j) . Simplify the above by making 
repeated use of the anti-symmetry of the epsilon tensor: 



! ffint(%A™(q) 



vsPur (q)) 



i(2n) 



(2g) 2 / dhdh S 3 (h + h)e^eP aXK f v {h)f a {l 2 ){i{\ - £))e^' 



+h")t 



(q + h)„{q + h) P e* SK (q)e s0 (q) {q - Z 2 ) M (g - l 2 ) p < {q)^M 



(q + h) 2 ~ qaqX (q + h) 2 + *e ' {q-hY ™* (q l 2 )* + ie _ 

l -^-{2g) 2 j dhdl 2 S 3 (h + h)e^ a0 e"^f u (h)f a (h)(i(l - ^)e^°+ l ^ t 



q a qx- 



X (^ur/3st;(q) ~ ^ S P^r(q)) 



(h)n(h) P z* SK (q)ts0(q) , (h)^h) P <p(q)t™(q) 



{q + hY qaq \q + hY+K ' (q~h) 2 ^(q-h) 2 +^e 



qaqx- 



i{2ir) 



(2g) 2 / dhdh 5 3 (h + l 2 )(i(l - e))(W~(q) " S vs Pur(q))e lih ° +h0)t 



x ( e^(h)^ u (h)q a ) ( eP° AK (h) P T a (h)qx 



0. 



eL(9)e s /3(?) 



e* sf3 (qYr K (q) 



(q + hY(q + h) 2 +ie (q - h) 2 - h) 2 + ie 



This expression vanishes because c'" tAk (l 2 ) p T a (h) = 0, as we have required in order to arrive at Eq. (|B-15|) . Therefore, 
the final evolution equations for the photon density matrix p will indeed be independent of the gauge parameter £. 
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